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THEORY  OF  THE  REACTION  SURFACE  WITH  Lef  *  1  AND  ITS 
APPLICATION  IN  THE  OXIDATION-SUBLIMATION  COUPLING  PROCESS 


Zhong  Jiakang 
(Zhejiang  University) 


i.  \  - 


ABSTRACT 


This  paper  presents  a  theoretical  analysis  on 
the  reaction  surface  of  a  locally  similar  laminar 
boundary  layer  with  addition  of  mass  and  chemical 
reaction.  Heat  transfer  and  mass  transfer  correla¬ 
tions  under  the  condition ^Lef  1  are  derived. 


As 


an  application  of  the  theory  of  the  reaction  surface 
the  oxidation  sublimation  coupling  process  on  the 
surface  material  is  investigated.  The  transition 
mechanism  and  transformation  condition  of  the  oxida¬ 
tion  and  sublimation  processes  from  rate-controlled 
to  diffusion-controlled  regimes  are  described.  The 
phenomenon  that  the  rate- sublimation  is  involved  in 
diffusion-controlled  oxidation  regime  is  shown.  Com¬ 
parisons  with  the  experimental  data  are  made  for  the 
case  of  graphite  stagnation  ablation,  and  quantita¬ 
tive  agreement  is  found.  (  W* '•  »'•  <  ?  <• 


SYMBOLS 

a_,  condensation  coefficient 

E 

C^  fractional  mass  of  a  unit  element 

Cp^  isobaric  specific  heat  of  a  unit  element 

C I  c.C  .  freezing  specific  heat 
pt  i  pi 

D.  equivalent  or  effective  two-dimensional  diffusion 

1  coefficient 

E  activation  energy 

f (n)  simplified  flow  function 

h^  £c^rfr  +  4j,  enthalpy  of  a  unit  element 

h?  formation  enthalpy  of  a  unit  element 

1  12 

h  h  +  -s-u  total  enthalpy 

s  • 


hTg  freezing  total  enthalpy 

k  surface  oxidation  rate  constant 

o 

S  y/pD.  Schmidt  number 

G  X 

Sc  !s,/(i  + 

T  temperature 

u,v,x,y  velocity  components  in  the  x,y  directions 


formation  rate  of  chemical  substance 
quantitative  relationship 

24ia 

dissociation  heat 
7y-l!k*ii+  reaction  heat 

h^-hg  dissolution  heat 

sublimational  fraction  of  a  unit  element  E 
viscous  coefficient 


v  -  -quantitative  coefficient 
p  mass  density 

Kf  conduction  coefficient  of  freezing  heat 

i  pu/ Peye 

X  1/(1  +  ctB) 


Lef  pDfCft/K,,  freezing  Lewis  number 

L£  hE~hE  •  sublimation  heat 

M.  molecular  weight  of  a  unit  element 

_i  r  ^ 

M  »•'  average  molecular  weight 


p  pressure 

Prf  freezing  Prandtl  number 

Prf  Prf/(1  +  aB) 

r  radius  of  revolution 

R  universal  gas  constant 

Rb  radius  of  the  head  of  the  stagnation  point 

E  chemical  element  of  the  sublimational  element  or 

surface  material 

e  condition  outside  the  fringe  of  the  boundary  layer 

eq  equilibrium  condition 

i  unit  element 

K  chemical  substance 

s  stagnation  point 


2 


'  derivative  against  n 

Superscript 

This  paper  received  on  March  17,  1982 

Many  references  (such  as  [1]  and  [2] ) concerning  the  theory 
of  the  reaction  surface  in  the  chemical  boundary  layer  have  been 
published.  In  prior  studies,  their  effort  mostly  considered 
the  characteristics  of  steady  flow  or  assumed  that  the  Lewis 
number-wars  1.  This  paper  uses  the  Blasius  function  to  solve  for 
the  characteristic  reaction  surface,  and  conducts  a  theoretical 
analysis  on  the  reaction  surface  of  a  locally  similar  laminar 
boundary  layer  with  addition  of  mass  and  chemical  reaction.  We 
derive  some  mass  transfer  and  heat  transfer  correlations  when 
Lef  *  1,  and  these  correlations  are  very  convenient  for  engineer¬ 
ing  applications.  As  an  application  of  the  theory  of  reaction 
surface,  we  conduct  a  study  into  the  coupled  oxidation-sublima¬ 
tion  process  on  the  surface  material  and  describe  the  transition 
mechanism  and  transformation  condition  of  the  oxidation  and  sub¬ 
limation  processes  from  rate  controlled  to  diffusion  controlled 
regimes.  The  phenomenon  that  the  rate  sublimation  is  included 
in  diffusion  controlled  oxidation  regime  is  shown.  Comparisons 
with  the  experimental  data  are  made  for  the  case  of  graphite 
stagnation  ablation,  and  quantitative  agreement  is  found. 

I .  Theory  of  the  reaction  surface 

We  introduce  the  Blasius  function  G  into  the  boundary  layer 
equation  and  use  G  as  the  new  independent  variable*  G  is  defined 


(1.1) 


Here 


g(,;  I.,  z>  -  «*(•»  I..  z)  5, 

G'(o;  z)  - 1  /jV'tt'*  in 


(1.2) 


assuming  that  the  gaseous  reaction  in  the  boundary  layer  is  con¬ 
fined  on  an  infinitely  thin  cross-section  (the  reaction  surface) 
and  the  flews  on  two  sides  of  the  reaction  surface  are  chemically 
frozen.  Based  on  this,  we  can  derive  the  following  reaction 
surface  equations  from  the  boundary  layer  equations  of  multi- 
component  reaction  gases  (3]  (j  *  0  is  the  horizontal  flow, 
j  =  1  is  the  axi-symmetrical  flow) : 

momentum  equation  f’  =  (1.3) 

diffusion  equation 

_  Cj  ****  Cim  _  _  St_  t*  (1.4) 

1  -  Gi  g;  Gf 

energy  equation 


Here 


Or,).-*?, 


(1.5) 


{ 


(1.6) 


(1.7) 

P'**\f*  p* 

and  G,  —  G(» ii  /.,  1),  G>  —  G(q;  S.)  and  G,  —  G(n\1mt  fr ,). 

There  are  existing  accurate  numerical  values  of  the  function 
G(n;  2).  that  we  can  readily  use  [4],  We  have  already  made 
the  following  assumptions:  (1)  the  flow  is  locally  similar; 

(2)  neglect  the  small  differences  of  the  specific  heat  between 
the  components;  (3)  L,  Pr^  and  Sc  are  all  constants;  (4)  neglect 
the  dissipation  term  in  the  energy  equation  (the  dissipitation 
effect  can  be  incorporated  in  the  computation  of  the  recovered 
enthalpy) ;  (5)  the  pressure  gradient  term  in  the  momentum 
equation  will  be  considered,  but  that  term  will  be  approximated 
as  0l(p,/p)  —  f*]  —  •Pit"  ;  here  a  is  a  selective  constant,  we  can 
set  a  =  1  +  2.6  (Tw/Tc)  [5]. 


By  using  the  exchanged  parameters  B  and  9^  of  the  intro¬ 
duced  mass  defined  as 


4 


u» »tiwt  m  wm w mw^rrwyrwTr mjrw;  rmr 


twt>vt 


’T™ST’VT»W'TJ! 


(1.8) 


(p»)»  “  «r<Oj  f„)B,  (p,*,).  “  6>(p9)m 

with  .  _ 

(1.9, 

We  can  write  the  mass  fraction  of  a  compound  on  the  surface  as 

(1.10) 


,  c?  +  Gj$jB 
...  .  .  J  +.GJB  . 


Combining  (1.4)  and  (1.8) ,  the  concentration  relation  on 
the  reaction  surface  is 

(1  +  GTBW  -M 1  -  G?)ftB  .  (1  -  c?)(l  ■+  OTBl  (1.Ha) 

"*  _  '  i  i  11 


1  +.B 


1  +  B 


cr 


Or  the  mass  fraction  of  the  chemical  element  is 


-  O  +  CTB)(c  r\  +  (1- 


Ck 


(1.11b) 


with 


1  +  B 


u  -  2  B*'f>  *  **  -  2 


The  expression  of  the  enthalpy  on  the  reaction  surface  can 
be  obtained  from  (1.4),  (1.5)  and  (1.10)  as 

At,  "  Gf(ir,)«  +  O  *”  C*)ir,  +  l*G* Ah^  (1  •  12) 

Here  we  have 

AA,»  —  2  A?[f„  —  (1  + 


/*-£r,Gr(l-G?)/Cf'(l-  C?) 


(1.13) 

(1.14) 


The  temperature  on  the  reaction  surface  can  be  computed 

* 

based  on  hT  .  This  is  to  say  that 

(1.15) 

By  using  (1.12),  the  heat  flux  from  the  boundary  layer  to 
the  surface  can  be  expressed  as 

-  qm  -  «4(0 if.)  [.(MX-  +  /*AA-  +  K0)B  2  A?  (r,„  -  B,)]  (1.16) 

Here  we  have 


B*(0;  ,.)  -  C(0;  P«,> 

P'i  V2| 


/(0)  -  Le,C’(0;  U,  l)/G’(0;fm,  Prf ) 


(1.17) 


(1.18) 


Based  on  the  definition  of  B  andthe  -f  boundary  condition 
correlation,  we  can  find  the  relationship  between  the  two  as 

w-/.  -  6\ 0 ;  1m±S')B/St  (1.19) 

During  the  computation  of  mass  transfer  and  heat  transfer, 
we  need  to  know  the  value  of  G'(0;/„  Z)'  .  For  engineering  appli¬ 
cations,  we  will  use  the  following  approximation  equation  to 
compute  the  injection  rate  (-f  )  and  its  influence  on  G  v®»  /«•»  ZT 
G'CO;  fm,  z )  -  G'( 0;  0»  2)  -  ^(-/.Z)  +  f,(-/.Z)*  d.20) 

Here 

...  G'(0;  o,  Z)  -  0A7Z*  (1.21) 

The  values  of  the  constants  and  are  listed  on  Table  1. 

They  are  related  to  Z.  After  detailed  numerical  comparisons 
[4],  this  equation  gives  very  high  accuracy.  So  the  substitution 
of  (1.19)_  into  (1.20)  can  give  a  second  order  equation  describing 
the  relationship  between  G'(°5/»*.z)  and  B  as 

«[  G'(0;  s,)V  +  AG'(0;  +  G’(o;  o,  5f)  -  0  (1.22) 

Here  a  =  £2B  and  b  =  “d+  CjB)  • 


TABLE  1. 


The  values  of  and  52 


>e 

• 

o 

N 

0.7 

6.8 

0.9 

1.0 

'  2.0 

.  e,  ...  0x770 

■  0.757 

0.753. 

■  0.745 

0.732 

0.702 

c,  ‘  ‘  0.0971  ' 

0. 163 

0.172 

0.T93 

0.192 

0.203 

G’(0;U,S ,) 

There  is  an  approximated  relationship  between^  G'(P»fi>  Ft,) 

.  G'(°;  tsjjsl.  —  LeJ  (1  23) 

6XQiL,Pr,)  ’  U*23) 

#  •- 

Here  the  exponent  n  is  related  to  the  injection  amount,  and 
their  relationship  is  given  by  the  table  below. 


*1) » \  nt '  i  "i".  \  i<«j  < 


Beside  this,  the  effect  of  density-viscosity  on  the  mass 
transfer  and  heat  transfer  is  expressed  by  +J~~i  —  \/pft/(tprfie')  . 

We  can  estimate  the  value  of  L  by  using  the  reference  temperature 
method,  or  to  simply  take  l  —.CfiwPml 

II .  Study  into  the  coupled  oxidation- sublimation  process 


By  using  results  from  the  above  study,  we  will  now  examine 
the  coupled  oxidation-sublimation  (or  evaporation)  process  of 
the  surface  material.  For  convenience,  we  will  specify  that  the 
oxygen  in  the  boundary  layer  will  only  form  EO  with  the  surface 
material  (expressed  as  chemical  element  E)  and  the  sublimated 
element  of  E  (may  be  vapor  of  a  single  molecule,  two  molecules 
or  multi -molecules)  will  only  form  two-molecule  EO  and  EN  with 
the  molecules  of  oxygen  and  nitrogen.  We  also  assume  that  the 
reaction  between  the  sublimated  element  and  air  is  a  fast  reaction 
(such  as  the  reaction  between  the  carbon  vapor  and  the  air) ; 
then  there  exists  a  complete  reaction  surface.  On  this  surface 
(fo  +  *o,)"“0  and  .  The  position  of  the  complete 

reaction  surface  is 

G*  -(* -  Se)/B(l  +  B.)  (2.1)  /31 

Here 

B .  -  (MeJMoXco  +  'o,\  +  (»/0  +  0>(tf*/W.0(Or  + 
and  (.Mcn/Mn)*  ™  <■£»/(<*  +  <■»,)*..  ,  *  is  the  carbon  consumption  of 

the  nitrogen- carbon  reaction. 

We  will  now  treat  the  oxidation  process  and  the  sublimation 
process  separately,  and  use  the  sublimation  ratio  A  (mass  loss 
due  to  sublimation/ total  mass  loss)  to  establish  the  linkage 
between  the  two.  Both  the  oxidation  and  the  sublimation  processes 
are  formed  by  two  basic  processes:  one  is  a  surface  chemical 
dynamical  process,  the  other  is  a  convection-diffusion  process 
in  which  either  the  reaction  material  (oxygen)  flows  toward  the 
surface  or  the  sublimated  material  flows  out  of  the  surface. 


o.-. 


. ,  v.y  v-v-v-v*. \*.y. 

JL  1  /Llk  \  .J\  1  ..  *.i  i_V.  V  .  V  -  V.  H 


Between  these  two  processes,  the  slower  one  of  the  two  will  be 
dominating.  Especially  when  the  rate  of  one  of  the  basic 
processes  is  far  less  than  that  of  the  other,  the  slower  one  of 
the  two  controls  the  entire  process  (oxidation  or  sublimation) . 
At  this  time,  the  rate  of  mass  loss  from  the  surface  actually 
equals  the  value  of  the  slower  basic  process. 


During  the  oxidation  process,  B  <'B#»  at  the  complete  reaction 
plane  on  the  surface.  The  rate  of  mass  loss  on  the  surface  can 
be  determined  by  a  combination  of  the  wall  dynamics  factor  and 
the  convection-diffusion  factor  within  the  boundary  layer.  It 
can  be  expressed  as 


V  Tfrrt  +  cn.\+(tli  +  cWt),1/0.±gL: 

■  ,0/flU)  . 


(2.2) 


Here  we  have 

_ a,*  —  (M£/M0)(l  +  By-%e-ua^UjMoyp\ 

<*,  -  (.Mt/McXl  -iXi  +  B)ae(0i  /.) 

They  are  the  coefficient  of  the  oxidation  rate  and  the  equivalent 
convection-diffusion  coefficient,  respectively.  The  former  char¬ 
acterize  the  speed  of  the  oxidation  reaction  of  the  surface 
material  under  a  certain  concentration  of  oxygen,  and  the  latter 
characterize  the  speed  of  transport  of  oxygen  to  the  surface 
which  is  a  complicated  function  related  to  the  boundary  layer 
factor. 


At  low  surface  temperatures  of  («a  «.«*,),  expression  (2.2)  can  lead 
to  the  limiting  case  of  oxidation  rate  control: 


£-*■0,  (c0  +  +  *<>,)»>  (pp)»  -*  (Mf- p  \  (c<,  +  coy  (2.3) 

'Mn  / 


8 


For  low  surface  temperature  (3t*  we  can  derive  the 

limiting  condition  for  the  diffusion-controlled  oxidation  as: 

(2.4) 

B  -+■  £•,  (c0  +  ‘ oX  — 0 >  • 

During  sublimation,  the  mass  loss  from  the  surface  as  deter 
mined  by  a  combination  of  the  surface  sublimation  dynamics 
factor  and  the  convection-diffusion  factor  is 

.1GMO--S777  wrwTl  C2 . 5) 

Here(«,)E  —  «eW»p/v 2xRMcTmis  called  the  coefficient  of  the  rate  of 
sublimation  for  the  sublimated  component  E.  It  represents  the 
speed  of  the  sublimation  process  on  the  surface  under  a  specific 
concentration  of  vapor; 


-  [(1  4-  C*)B/C* ]a,(0;  O  -  [(1  +  B)B/(B  -  B.)]«,( 0;  /„) 

is  an  equivalent  convection-diffusion  rate  coefficient;  it  repre¬ 
sents  the  speed  of  disappearance  of  the  sublimated  material  from 
the  surface  and  it  is  also  a  complicated  function  of  the  boundary 
layer  factor. 


When  the  surface  temperature  is  not  too  high,  we  have  ((a,)E  « <**,) 
B —*  Bt,  (i.m  —*■  0,  (2.6) 

This  is  a  limiting  condition  for  the "rate-controlled  sublimation. 

In  the  region,  the  complete  reaction  plane  is  still  maintained  on 
the  surface;  the  sublimated  element  is  instantaneously  burned 
completely  on  the  surface.  The  oxygen  that  is  transported  to 
the  surface  is  partially  consumed  by  the  sublimated  element,  and 
the  remaining  part  then  forms  the  oxidation  reaction  with  the  solid 
wall.  The  mass  loss  on  the  surface  equals  the  value  of  diffusion- 
controlled  oxidation.  This  indicates  that  during  the  coupled 
oxidation-sublimation  process,  the  rate  controlled  sublimation 
is  contained  in  the  region  of  diffusion  controlled  oxidation, 
and  the  rate  controlled  part  of  the  sublimation  increases  with 
increasing  temperature,  until  the  surface  mass  loss  is  completely 
due  to  sublimation  (X  =  1) .  The  sublimation  ratio  can  be  deter¬ 
mined  by  the  equation  below 

i  -  £  /-)*.  (2.7) 


When  the  surface  temperature  is  sufficiently  high  (( <*,)«»«/,) 
we  have 


(pv)m  tm)B 


(2.8) 


This  is  the  diffusion-controlled  sublimation  regime.  This  is 
equivalent  to  the  condition  that  the  wall  surface  maintains  a 
limiting  equilibrium  sublimation  condition  when  the  transport 
in  the  boundary  layer  is  very  slow  so  that  the  sublimated  gases 
cannot  flow  away  from  the  surface  fast  enough.  At  this  time, 


the  solid  wall  oxidation  reaction  no  longer  exists  (X  =  1)  since 
the  complete  reaction  plane  has  gone  away  from  the  surface. 

Ill .  Results  and  conclusions 

We  will  now  discuss  the  oxidation  and  sublimation  of  carbon 
material  as  an  example.  Figure  1  shows  the  variation  of  the  mass 
exchange  parameter  B  with  surface  temperature  Tw  during  a  non¬ 
equilibrium  oxidation  process.  The  surface  oxidation  dynamics 

data  are  obtained  for  (1)  "slow  reaction":  n  =  i,  E  =  42.3 

42  ^  1/2* 

Kcal/g  molecule,  k  =  2.18  x  10  g/cm  S  (atmospheric  pressure)  ' 

°  1 

(2)  "fast  reaction":  n  =  •=■,  E  =  44.0  Kcal/g  molecule,  k  =3.28 

82  ^  1/2  ^ 

x  10  g/cm  S  (atmospheric  pressure)  '  .  The  former  corresponds 

to  heat  dissolution  graphite,  and  the  latter  corresponds  to  common 
industrial  graphite  [7].  By  defining  a  Damkoheler  number  = 

=  ad.j/ach,  the  transformation  condition  of  the  controlled  regime  can 
be  determined.  For  example,  we  can  select  Dm  =  100  and  0.01 
(corresponding  to  B  values  of  0.00173  and  0.171,  respectively)  as 
the  boundary  values  for  the  non-equilibrium  transition  regime. 

For  the  flow  around  the  stagnation  point  under  the  condition  that 
Rjj  =  10  Cm,  the  corresponding  temperature  boundary  values  for  the 
"fast  reaction"  are  770K  and  1000K,  respectively.  For  "slow 
reaction"  the  corresponding  values  are  HOOK  and  1700K  respectively 
(Figure  1) .  Under  the  special  condition  of  flow  around  the  stag¬ 
nation  point,  the  boundary  value  of  temperature  is  determined  by 
the  radius  of  the  head  and  the  oxidation  dynamics  data.  Since  n  = 
1/2,  pressure  has  nothing  to  do  with  it.  We  can  select  the  air 
concentration  outside  the  fringe  of  the  boundary  layer  to  be  (cq  + 
c  )  =  0.23  and  (c„  +  c„  )  =  0.77  in  our  computation.  Figure  2 

shows  the  B  versus  T  curves  for  the  equilibrium  as  well  as  the 
non-equilibrium  sublimation  process.  The  major  compositions  of  the 
carbon  vapor  are  assumed  to  be  C,  C2,  C3,  C4  and  C 5«  The  equili¬ 
brium  concentration  of  the  sublimated  element  can  be  compiled  from 
the  equilibrium  vapor  pressure  pE  as  cE  =  M^jPg  e^/Mp.  The 
equilibrium  vapor  pressure  of  each  carbon  element  is  also  given 


m  the  JANAF  thermodynamics  table  [8],  and  can  be  computed  by  the 
following  equation 

log  MpZxf  ”  «  +  ^/7*(K)  £2 

The  constants  a  and  b  are  listed  in  Table  2.  Table  2  also  shows 
the  condensation  coefficients  of  the  carbon  vapor  aE.  These 
coefficients  are  obtained  from  [9]. 

Figure  2  shows  a  curve  with  the  assumption  that  a_  =  1.  We 

£i 

can  see  that  the  aE  =  1  results  is  very  similar  to  the  equilibrium 
value,  while  the  deviation  is  larger  when  a£  is  obtained  from  the 
values  in  Table  2. 


TABLE  2.  The  condensation  coefficient  of  the  carbon  vapor  aE 
and  the  constants  a  and  b  in  the  equilibrium 
constant  equation 


E  - 

—  c 

c. 

c. 

c. 

c. 

«« 

.  0.24 

o.so 

0.023 

0.25 

0.0019 

4 

8.096 

9.607 

9.627 

10.013 

10.464 

'  * 

-3.709 

-4.262 

-3.960 

-4.889 

-4.927 

(§*¥***(*»- io 

•s  ■  1 - 

**2* - 


Figure  1. 
B  with  T_ 


w 


The  variation  of 
in  an  oxidation 


process , 

1 —  mass  exchange  parameter  B; 

2—  diffusion-controlled  oxida¬ 
tion  regime;  3 — "fast  reaction" 
(r.=10Cm)  ;  4 — "slow  reaction" 
(R?*10Cm) ;  5 — surface  temper¬ 
ature  T  °K 

w 


Figure  2.  The  variation  of 

B  with  T  in  a  sublimation 
w 

process . 

1 —  mass  exchange  parameter  B; 

2 —  p=0,k  atmospheric  pressure; 

3—  diffusion-controlled  oxida¬ 
tion  regime;  4 — surface  temp.; 

5 — equilibrium  sublimation;  6— 
non-equilibrium  sublimation; 

7 — a_  obtained  from  values  i-n 
Tabll  2 


Based  on  the  energy  equilibrium  relation  at  the  gas-solid 
interface  and  (1.16) ,  the  heat  flux  into  the  interior  of  the  . 
solid  can  be  expressed  as 

q »  —  «*(0;  /.)[(*r,)»  -  +  lm*Ja  -  I(0)BL ,]  +  aq,t  -  «rTl  (3.2) 


Here  we  have 

£ji,t  -  -  E  Ue,,X  -  (1  +  +  ‘O' 

Cy 

—  (1  +  B)(cj  -h  cjt)m]&Qcj  +  E  K1  +  b)ce.w  “ 

I*C| 

Z,-(l  ”  l)L,  — 


*-c 


(3.3) 


(3.4) 


TABLE  3.  Reaction  heat  (cal/g) 


®  nuKSft 

a0o,  -  3700 

a0c,  —  6100 

tc«H200  2c,«460Q 

“  »000 

a  0c,  -  9200 

£c,-81Q0  £c,-3700 

a  Qco  *■  1600 

a  0c,  -  9600 

‘  * 

a0c5>  13500 

A0C, -.10500  . 

-  5000 

1 — gas  phase  reaction  heat;  2 — sublimation  heat 


/33 


The  quantities  aq„  and  ^T*,  represent  the  energy  loss  /34 

due  to  the  energy  transmitted  to  the  wall  by  the  heated  gas 

radiation  and  the  outward  radiation  from  the  surface  of  the  wall, 

respectively.  The  various  types  of  reaction  heat  are  listed  in 

Table  3.  From  (3.2)  we  can  see  that  the  reaction  plane  can 

* 

directly  influence  the  value  of  q^  through  the  factor  I  Ah^ 

When  I  =  1(0),  or  when  the  reaction  plane  is  in  the  wall  surface, 
(3.2)  can  be  simplified  into  the  result  of  Lees  [6]  for  the 
frozen  boundary  layer. 


Finally,  in  order  to  test  the  accuracy  of  our  theory,  we 
have  conducted  computation  on  the  ablation  of  graphite  at  the 
axial-symmetric  stagnation  point  in  dissoluted  air.  The  heat 
transfer  coefficient  of  the  dissoluted  air  a^(0.0)  for  the  flow 
around  the. stagnation  point  can  be  approximated  as 

«*(0;  0)  —  0.113  V p,l  g/an^-s  (3.5) 
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Here  pg  is  in  atmospheric  pressure,  is  in  Cm.  Due  to  the 
unstable  nature  of  the  transport  coefficient,  we  select  the 
following  two  values  —  0.7  and  ^  (•$*)»,  “  ;  (£)  S,  —  0.965 

and  Le^  =  1,  and  assume  that  the  sublimation  is  in  equilibrium. 

The  comparison  of  the  theoretical  values  and  experimental  values 
is  shown  in  Table  4.  Under alcw  pressure  environment,  the  graphite 
is  basically  in  thermal  chemical  ablation,  and  we  can  see  that 
there  is  a  quantitative  agreement  between  the  theoretical  value 
and  the  experimental  value. 

By  using  the  turbulent  transport  coefficient  to  replace  the 
laminar  transport  coefficient  in  the  boundary  layer  equation  and 
use  the  Reynolds  approximation,  this  method  can  be  expanded  to 
the  turbulent  boundary  layer. 

TABLE  4.  Comparison  between  theoretical  value  and  experimental 
value . 


9n 

aim. 

T. 

K 

<p  <P-).v/£7ir  (S/**'XM*/*W'* 

Q>  m  «  a 

<•)  • 

(b)  - 

0.12 

2300 

0.0140 

0.0174 

0.017C2520-2620K) 

2600 

0.0141 

0.0175 

0.12- 

3600 

0.0731 

0.0931 

0.069-0.077  (3650K) 

3700 

0.0103 

0.0993 

0.122-0.135  (3700K) 

0.32  ' 

2300 

0.0140 

0.0173 

0.01S  (2520— 2700K) 

2700 

0.0141 

0.0175 

0.30 

3000 

0.0130 

0.0115 

0.017  (3000 — 3200K) 

3200 

.  0.0169 

0.0209 

1 

0.32 

3600 

0.03(9 

0.0412 

0.037-0.052  (3600K) 

3700 

0.0373 

0.0709 

0.062-0.066  (3710K) 

CO  l^*<*0*,  “  0.7  ft  5c<fc(Sc)w,  •  0.965 
(b)  ie-  0.965 


2  1/2 

i — (g/cm  )  (Cm/atmospheric  pressure)  ;  2— theoretical  value; 
3 — measured  value 


The  first  draft  of  this  paper  was  written  in  1976.  The  author 
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JTHEORY  OF  THE  REACTION  SURFACE  WITH  Let\\ 

^ AND  ITS  APPLICATION  TO  ^  ^ 

OXIDATION-SUBLIMATION  COUPLING  PROCESS 


Abstract, 


This  paper  presents  a  theoretical  analysis  of  the  reaction  surface  for  the  locally 
similar  laminar  boundary  layer  with  addition  of-  mass  and  chemical  reactions.  Beat- 
transfer  and  mass-transfer  correlations  for  Lct  *v  1  are  obtained.  As  an  application 
the  oxidation-sublimation  coupling  process  is  investigated.  The  transition  mechanism 
from  rate-controlled  to  ' diffusion-controlled  regimes  is  described.  The  Phenomenon  that 
.the  rate-sublimation  is  involved  in  diffusion-controlled  oxidation  regime  is  displayed. 
Comparisons  with  the  experimental  data  are  made  for  the  ease  of  graphite  stagnation 
ablation  and  quantitative  agreement  is  found. 
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